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ACCELERATION AND SUBSTRUCTURE CONSTRAINTS IN A QUASAR OUTFLOW 
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ABSTRACT 

We present observations of probable line-of-sight acceleration of a broad absorption trough of C IV 
in the quasar SDSS J024221. 87+004912. 6. We also discuss how the velocity overlap of two other 
outflowing systems in the same object constrains the properties of the outflows. The Si IV doublet 
in each system has one unblended transition and one transition which overlaps with absorption from 
the other system. The residual flux in the overlapping trough is well fit by the product of the residual 
fluxes in the unblended troughs. For these optically thick systems to yield such a result, at least one of 
them must consist of individual subunits rather than being a single structure with velocity-dependent 
coverage of the source. If these subunits are identical, opaque, spherical clouds, we estimate the cloud 
radius to be r ~ 3.9 x 10^^ cm. If they are identical, opaque, linear filaments, we estimate their 
width to be w ~ 6.5 x 10^^ cm. These subunits are observed to cover the Mg II broad emission line 
region of the quasar, at which distance from the black hole the above filament width is equal to the 
predicted scale height of the outer atmosphere of a thin accretion disk. Insofar as that scale height is 
a natural size scale for structures originating in an accretion disk, these observations are evidence that 
the accretion disk can be a source of quasar absorption systems. Based on data from ESO program 
075.B-0190(A). 

Subject headings: quasars: general, absorption lines, individual (SDSS J024221.87+004912.6) 



1. INTRODUCTION 

Absorption systems in active galactic nuclei (AGN) can 
be classified as intrinsic, which are associated with the 
active nucleus and often outflowing from it, and interven- 
ing, which originate from clouds external to the AGN's 
environment. Determining an absorption system's clas- 
sification can be difficult. A reliable indication of intrin- 
sic absorption is time variability, such as a shift in the 
velocity of a given feature or changes in its absorption 
strength as a function of velocity. In the broad absorp- 
tion line (BAL) troughs > 1000 km s~^ wide which are 
most often found in the luminous AGN known as quasars, 
reports of time variable absorption strengths have been 
relatively common (Bromage et al. 1985; Foltz et al. 
1987; Voit, ShuU, & Begelman 1987; Smith & Penston 
1988; Turnshek et al. 1988; Barlow, Junkkarinen, & 
Burbidge 1989, 1992a; Barlow et al. 1992b; Barlow 1994; 
Hamann et al. 1995; Michalitsianos, Oliversen, & Nichols 
1996; Hamann, Barlow, & Junkkarinen 1997b; Barlow, 
Hamann, & Sargent 1997; Hamann et al. 1997a; Hall 
et al. 2002; Gallagher et al. 2004; Misawa et al. 2005; 
Lundgren et al. 2006). Such variability can even include 
the appearance or disappearance of absorption systems 
(Koratkar et al. 1996; Ganguly, Charlton, & Eracleous 
2001; Ma 2002; Gallagher et al. 2005; Leighly et al. 
2005). In contrast, velocity shifts in BAL outflows have 
been reported in only Q 1303-1-308 (Vilkoviskij & Irwin 
2001) and Mrk 231 (Rupke, Veilleux, & Sanders 2002, 
and references therein), although Gabel et al. (2003) 
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have observed deceleration of a narrow absorber in the 
Seyfert 1 NGC 3783. 

Acceleration must occur for AGN outflows to reach 
their observed velocities. Nonetheless, velocity shifts in 
AGN outflows are seen quite rarely because acceleration 
of an AGN outflow does not automatically translate into 
a change in its observed velocity profile, and vice versa. 
For example, a fixed mass loss rate into an outflow with a 
time-invariant driving force would yield a time-invariant 
acceleration profile with distance in the outflow, and 
thus produce unchanging absorption troughs. Arav ct al. 
(1999) illustrate how radial acceleration of gas crossing 
our line of sight with a non-negligible transverse velocity 
produces an observed absorption trough with a broad- 
ened radial velocity profile that does not change with 
time. Since our lines of sight to AGN are essentially ra- 
dial, and since AGN are fed by accretion disks consisting 
of gas with predominantly orbital velocities, most AGN 
outfiows are expected to have non-negligible transverse 
as well as radial velocities. Thus, most intrinsic absorbers 
likely are exhibiting acceleration, disguised as a trough 
broader than the thermal or turbulent velocity width of 
the gas. 

What are we then to make of cases where an outflow 
does exhibit a velocity shift? First, note that when our 
line of sight intersects the origin of an outflow, the ab- 
sorption trough can start at zero line-of-sight velocity in 
the AGN rest frame, at least for ions present at the origin 
of the outflow. Ions present only downstream in an out- 
flow, or lines of sight intersecting an outflow only down- 
stream from its origin due to curvature in the flow lines, 
will produce 'detached' absorption troughs which do not 
start at zero velocity, as will a shell of material ejected 
in an intermittent outflow. With that in mind, consider 
possible explanations for a velocity shift observed in a 
detached absorption trough. Such a shift can be pro- 
duced by changes in the ionization state as a function of 
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velocity in a fixed outflow, by changes in the acceleration 
profile or geometry (or both) of such an outflow due to 
changes in the driving force or mass loss rate, or by ac- 
tual line-of-sight acceleration of a shell of material from 
an intermittent outflow. Observations of velocity shifts 
are therefore worthwhile because they may yield insights 
into specific scenarios for quasar absorbers. 

Here we present multiple-epoch observations (§2) of a 
quasar in which a broad absorption line trough of C IV 
increased in outflow velocity over 1.4 rest-frame years 
(§3). We also discuss how two overlapping outflows in 
the same quasar provide constraints on the properties of 
those outflows (§4). We end with our conclusions in §5. 

2. OBSERVATIONS 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) 
is using a drift-scanning camera (Gunn et al. 1998) on a 
2.5-m telescope (Gunn et al. 2006) to image lO'^deg^ of 
sky on the SDSS ugriz AB magnitude system (Fukugita 
et al. 1996; Hogg et al. 2001; Smith et al. 2002; Pier 
et al. 2003; Ivezic et al. 2004). Two multi-fiber, double 
spectrographs are being used to obtain resolution R ~ 
1850 spectra covering ~3800-9200 A for -10^ galaxies to 
r = 17.8 and ~10^ quasars to i = 19.1 (i = 20.2 for z > 3 
candidates; Richards et al. 2002). 

The Zem = 2.062 BAL quasar SDSS 
J024221.87-)-004912.6 (Schneider et al. 2002; Re- 
ichard et al. 2003; Schneider et al. 2005; Trump et al. 
2006), hereafter referred to as SDSS J0242+0049, was 
observed spectroscopically three times by the SDSS (Ta- 
ble [T]). We selected it for high- resolution spectroscopic 
foUowup because of the possible presence of narrow 
absorption in excited-state Si II and C II at z = 2.042. 
A spectrum obtained with the ESO Very Large Tele- 
scope (VLT) Unit 2 (Kueyen) and Ultra- Violet Echelle 
Spectrograph (UVES; Dekker et al. 2000) confirms the 
presence of narrow, low-ionization absorption at that 
redshift,^ analysis of which will be reported elsewhere. 

We observed SDSS J0242-h0049 with UVES on the 
VLT UT2 on the nights of 4-5 September 2005 through a 
1" slit with 2x2 binning of the CCD, yielding R ~ 40000. 
The weather ranged from clear to thin cirrus, with 
0.8 - 1.0" seeing. SDSS J0242-I-0049 was observed for a 
total of 5.75 hours in two different spectral settings, yield- 
ing coverage from 3291-7521 A and 7665-9300 A. Each 
exposure was reduced individually with optimum extrac- 
tion (Horne 1986), including simultaneous background 
and sky subtraction. Telluric absorption lines were re- 
moved for the red settings using observations of telluric 
standard stars. A weighted co-addition of the three ex- 
posures of each spectral setting was performed with re- 
jection of cosmic rays and known CCD artifacts. Finally, 
all settings were rebinned to a vacuum heliocentric wave- 
length scale, scaled in intensity by their overlap regions, 
and merged into a single spectrum with a constant wave- 
length interval of 0.08 A (Figure[l]). The SDSS spectra aU 
share a common wavelength system with pixels equally 
spaced in velocity, and so for ease of comparison we cre- 
ated a version of the UVES spectrum binned to the those 
same wavelengths but not smoothed to the SDSS resolu- 
tion. 

* The weak, narrow character of that absorption led to the classi- 
fication of this object as a high-ionization BAL quasar by Reichard 
et al. (2003) and Trump et al. (2006) based on its SDSS spectrum. 



3. BROAD ABSORPTION LINE TROUGH VELOCITY 
SHIFTS 

The broadest absorption lines in SDSS J0242-f-0049 oc- 
cur at a redshift z ~ 1.87988 {v = -18400 km s^^ relative 
to the quasar) in Lya, N V, Si IV and C IV (Figure [l]). 
There is an offset between the peak absorption in C IV 
and Si IV. The redshift z = 1.87988 was determined from 
the deepest absorption in the Si IV trough, and does not 
match the deepest C IV absorption. This can be ascribed 
to a changing ionization state in the outflow as a function 
of velocity. 

Comparison of the SDSS and UVES spectra suggested 
a shift in the position of the C IV trough at this redshift. 
To investigate further, continuum regions around that 
trough and the Si IV trough at the same redshift were 
fitted and used to normalize all observed spectra. (The 
Lya and N V troughs lie outside the SDSS wavelength 
range.) For each epoch, the C IV and Si IV regions were 
fit separately with third order Legendre functions using 
SPLOT in IRAF.^ The continuum sample windows were 
selected to avoid emission lines in the quasar rest frame 
(Vanden Berk et al. 2001). 

The extent of any shift can be measured by minimizing 
the between the normalized pixel-by-pixel fluxes in 
the spectra when shifted by an integer number of pixels 
m (assuming pixels equally spaced in velocity): 

2 _ 1 (/2,t - /l/t+m)^ ..N 

where /2,i and /i,i+m represent the flux in spectra from 
epochs 1 and 2 at pixels i and i -I- m, respectively, N is 
the total number of pixels extracted from each spectrum 
for comparison and a is the error for the flux at each 
pixel. 

The SDSS spectra from epochs 51821 and 52188*^ were 
compared with the UVES spectrum from epoch 53619 
(Table [1]) . A clear shift was found in C IV and a po- 
tentially smaller shift in Si IV. Neither trough shows a 
detectable shift between the SDSS spectra from epoch 
51821 and epoch 52188, and neither would be expected 
to do so if the observed long-term shift was due to a con- 
stant acceleration (the shift between those two epochs 
would be < 0.5 pixel for C IV). In light of this, the 
test was conducted again, using a weighted average of 
all three SDSS spectra, with mean epoch 52066. From 
that comparison we conclude that the shift in C IV is 
3 ± 1 pixels with 95.4% confidence {2a). Zero velocity 
shift in C IV can be excluded with 99.9998% confidence. 
For Si IV, the shift is 1 ± 3 pixels at 95.4% confidence. 
Plots of these spectra are shown in the top two panels of 
Figure [3l It is important to note that there is no shift in 
the nearby narrow absorption lines. Also, both troughs 
appear to keep a relatively constant intensity, within the 
uncertainties. The bottom panel of Figure [3] shows the 
excellent match to the epoch 53619 UVES spectrum that 

^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatories, which 
is operated by AURA, Inc., under contract to the National Science 
Foundation. 

Since the SDSS spectra from MJD 52177 and MJD 52199 are 
noisier at short wavelengths than the SDSS spectrum from MJD 
51821 and since visual inspection of them revealed no obvious dif- 
ference in their BAL troughs, a weighted co-add of their spectra 
was made, with mean epoch 52188. 
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results when the epoch 52066 average SDSS spectrum is 
shifted by 3 pixels. 

The middle panel of Figure [3] may suggest that the 
long-wavelength end of the C IV trough has a greater 
shift than the short-wavelength end. Splitting the C IV 
trough into two sections, we find that is minimized at 
a shift of 2^1 pixels for the short-wavelength end and a 

shift of pixels for the long-wavelength edge, but that 
a uniform shift produces a marginally lower minimum 
overall x^. Thus, while there is possible evidence for a 
nonuniform velocity shift of the C IV BAL trough, the 
current data are of insufficient quality to prove its exis- 
tence. Many physical effects could produce a nonuniform 
shift (expansion of an overpressured, accelerated shell of 
gas from an intermittent outflow, to give one example). 

A shift of one SDSS pixel corresponds to a velocity shift 
of 69 kms^^ in the observed frame or 22.5 kms~^ in the 
quasar rest frame (z = 2.062). A shift of 3 ± 1 SDSS 
pixels (2ct) over a rest-frame time span of 1.39 years thus 
gives an acceleration of a = 0.154 ± 0.025 cm s~^, where 
the error is la. Previously claimed accelerations for BAL 
troughs are much lower than that, at a = 0.035 ± 0.016 
cm s~^ over 5.5 rest-frame years in Q 1303-1-308 (Vilko- 
viskij & Irwin 2001) and a = 0.08 ± 0.03 cm s"^ over 
12 rest-frame years for Mrk 231 (Rupke et al. 2002). 
Our observation is more similar to that of Gabel et al. 
(2003), who determined the deceleration of C IV, N V 
and Si IV in a narrow absorption system in a Seyfert 
galaxy and found (for C IV) relatively large values of 
a = -0.25 ± 0.05 cm s'^ and a = -0.10 ± 0.03 cm s'^ 
over 0.75 and 1.1 rest-frame years, respectively. All of 
those observations involved much narrower troughs than 
is the case in SDSS J0242-I-0049. Also, the la rela- 
tive uncertainty associated with the acceleration of SDSS 
J0242-f0049 is lower than the previous BAL measure- 
ments. These factors make SDSS J0242-H0049 a robust 
case for line-of-sight acceleration of a true BAL trough. 
Still, it should be kept in mind that all these accelerations 
are much smaller than the a ~ 100 cm s~^ predicted for 
the main acceleration phase of a disk wind in the model 
of Murray et al. (1995). 

Furthermore, BAL troughs can vary for several rea- 
sons. These include acceleration or deceleration along 
the line of sight of some or all of the absorbing g as , a 
change in the ionization state of some or all of the gas, 
or a change in C{v) — the covering factor of the gas as a 
function of the line-of-sight velocity — due to the move- 
ment of gas into or out of our line of sight, for example 
due to a change in flow geometry (see the introduction 
and §3.3 of Gabel et al. 2003). In many cases of vari- 
ability all of the above origins are possible, but there 
are cases where acceleration is very unlikely to be the 
cause (see below). Because of this, to be conservative 
we cannot assume that BAL trough variability is due to 
acceleration even though acceleration could be the cause 
of much of the observed variability. 

Fig. 2 of Barlow et al. (1989) and Fig. 2 of Barlow 
et al. (1992b) are cases where observed time variabil- 
ity of BAL troughs is almost certainly due to a change 
in the column densities of an ion at certain velocities 
(whether due to a changing ionization or to bulk motion 
into the line of sight), not due to a given ionic column 
density changing its velocity. More ambiguous cases are 



illustrated by C IV in Q 1246-057 (Fig. 3 of Smith & 
Penston 1988) and Si IV in Q 1413-f 117 (Fig. 15 of Turn- 
shek et al. 1988). In both of those cases, a second-epoch 
spectrum shows more absorption at the short-wavelength 
edge of the trough in question. That could be because 
gas at lower outflow velocities in the trough was acceler- 
ated to higher velocities. Yet in both cases, the trough 
away from the short-wavelength edge is unchanged be- 
tween the two epochs. If acceleration was the cause of 
the variability, a reduction in covering factor or optical 
depth, or both, might be expected at the lower velocities 
where the gas originated. No reduction is seen, arguing 
against the line-of-sight acceleration hypothesis for these 
cases of trough variability. 

While every case for acceleration in a BAL trough will 
be ambiguous at some level, comparing the variability 
we report in SDSS J0242-f0049 to previous cases leads 
us to believe that ours is the least ambiguous case seen to 
date of acceleration in a true BAL trough (> 1000 kms~^ 
wide). Monitoring the future behavior of the z = 1.87988 
absorption in this quasar would be very worthwhile, to 
see if the acceleration was temporary, is constant, in- 
creasing, or decreasing, or varies stochastically. The lat- 
ter might occur if the velocity shift is due to a variable 
flow geometry or to ionization variations as a function of 
velocity caused by a fluctuating ionizing luminosity. (Re- 
call from Figure 2 that this system shows some evidence 
for ionization stratification with velocity, in the form of 
an offset between the velocities of the peak Si IV and C IV 
absorption.) As this quasar is located in the equatorial 
stripe of the SDSS, which has been repeatedly imaged 
over the past 7 years, it should eventually be possible to 
search for a correlation between its ultraviolet luminosity 
and the acceleration of this system. (From the spectra 
alone, there appears to be a 5-10% increase in the lumi- 
nosity of the object over the time spanned by the three 
SDSS spectra, but no information is available on longer 
timescales since the UVES spectrum is not spectropho- 
tometrically calibrated.) BAL trough velocity shifts are 
also expected if BAL quasars are a short-lived phase dur- 
ing which material is expelled from the nuclear region 
(Voit, Weymann, & Korista 1993). In such a model the 
accelerating trough in SDSS J0242-I-0049 could be inter- 
preted as gas unusually close to the quasar, currently 
experiencing an unusually large radiative acceleration. 

4. OVERLAPPING SI IV TROUGHS 

There is a possible case of line-locking involving Si IV 
in SDSS J0242-H0049. Stable fine-locking in a given dou- 
blet occurs wfien two conditions are met. First, the ve- 
locity separation between two absorption systems at dif- 
ferent redshifts must be very nearly equal to the velocity 
separation of the two lines of a doublet seen in both sys- 
tems (Braun & Milgrom 1989). Second, the reduction in 
line-driven acceleration of tfie shadowed system due to 
the reduced incident fiux in one component of the dou- 
blet must result in its acceleration being the same as that 
of the shadowing system. This latter condition may be 
difficult to meet in AGN outflows, where many lines con- 
tribute to the radiative acceleration and there may also 
be substantial non-radiative acceleration. Nonetheless, 
some spectacular examples of apparent line-locking in 
AGN do suggest that it can in fact occur (e.g., Srianand 
et al. 2002), even if only rarely. 
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As shown in Figure H in SDSS J0242+0049 there is 
narrow Si IV absorption a,t z — 2.0476 (hereafter system 
A') and a broad Si IV trough centered at about z — 2.042 
(hereafter system A). Si IV line- locking of a third absorp- 
tion system to system A' or A would result in absorption 
1931 kms~^ shortward of those redshifts, at z = 2.0280 
or z = 2.02245 respectively. What is observed in the 
spectrum, however, is broad absorption in between the 
expected redshifts, centered at z = 2.0254 (hereafter sys- 
tem B). Both systems are observed in other transitions 
as well, with system B having more absorption in N V 
and C IV but less in S IV and Mg II. 

In this section we consider first the optical depths and 
covering factors of these overlapping systems, with in- 
triguing results. We then consider whether they could 
be line-locked or in the process of becoming line-locked. 

4.1. Si IV Trough Optical Depths and Covering Factors 

It is useful to determine if the Si IV troughs un- 
der consideration are optically thick or not. Figure [5] 
shows the absorption profiles in velocity space relative 
to z ~ 2.0476 or to the corresponding line-locked red- 
shift of z = 2.0280. System A+A', seen unblended in the 
bottom panel, is free from contamination in the blended 
trough (middle panel) at —900 < v < —650 kms^^. At 
those velocities, absorption from the A1402 component 
of System A-|-A' (bottom panel) appears so similar in 
shape and intensity to absorption from the intrinsically 
stronger A1393 component (middle panel) that we can 
conclude system A+A' is optically thick in Si IV. For 
system B (seen unblended in the top panel) we must see 
how well various combinations of optical depth, cover- 
ing factor, and geometry (Rupke, Veilleux, & Sanders 
2005) can reproduce the profile of the trough composed 
of blended absorption from system B and the optically 
thick system A-l-A' (middle panel). 

For an unblended doublet, at each velocity v the nor- 
malized residual intensities Ii and I2 (in the stronger and 
weaker lines, respectively) can be related to the optical 
depth in the stronger transition r and the fraction of the 
emitting source covered by the absorber along our line of 
sight, the covering factor C (e.g.. Hall et al. 2003): 

/i(i;) = l-a(l-e-^") (2) 

J2(«)-l-a(l-e-^"") (3) 

where ii measures the relative optical depths of the lines. 
For the Si IV AA1393, 1402 doublet, R = 0.5. In each ab- 
sorption system we have only one unblended component, 
but it can still be used to model the other component. 
(For comparison, the two unblended troughs are over- 
plotted on the blended trough in the top panel of Fig- 
ure [6l) 

First we test whether system B can be optically thin, 
with Ct, = 1. Using this assumption and equations [2] and 
[3l the optical depth Tu(A1402, B) was calculated from the 
observed trough of Si IV A1393 in system B. The blended 
trough profile in this model should be exp[— t„(A1402, B)] 
times the profile of Si IV A1393 in system A-^A'. (The 
latter profile is taken as identical to the A1402 trough pro- 
file at z = 2.0476 since system A+A' is optically thick.) 
The resulting model blended-trough profile is compared 
to the observed blended-trough profile in the second 



panel of Figure [6l Optically thin absorption from sys- 
tem B falls short of explaining the depth of the blended 
trough. 

Next we test whether system B can be extremely op- 
tically thick, so that the depth of its absorption is deter- 
mined only by Cy. In this case, we have two absorption 
systems absorbing at each v, but with different C^. The 
total absorption is determined by C^^uended, which de- 
pends on what parts of the emitting source (s) are covered 
by neither absorption system, by just one, or by both. 
That is, the total absorption depends on the extent to 
which the two systems overlap transverse to our line of 
sight and cover the same parts of the source. We can rule 
out the limit of minimum overlap, which yields maximum 
coverage of the source: C^^uended = min(Cyi + Cb, 1). In 
that case Ca + Cs > 1 at all w, but we do not observed 
Cv,biended = 1 at all V. Another limiting case is maxi- 
mum overlap of the absorption systems, which minimu- 
mizes the source coverage: Cv,biended — max(CA,CB). 
The results of that model are shown in the third panel 
of Figure [51 It is not an improvement over the opti- 
cally thin model. However, at almost all velocities the 
maximum-overlap model has more residual flux than seen 
in the data, while the minimum-overlap model has less. 
Thus, overlap in Cy which is less than the maximum 
possible by a velocity-dependent amount can explain the 
data. Such spatially-distinct, velocity-dependent partial 
covering has been seen before in other quasars (see the 
Appendix to Hall et al. 2003). 

The last case we consider is one where each covering 
fraction describes the fractional coverage of the other ab- 
sorption system as well as of the continuum source, so 

that Iv,blended = IaIb and Cv^Mended = Ca + Cb—CaCb 

(this is case 3 of Rupke et al. 2005). The results of this 
model are shown in the bottom panel of Figure [51 again 
assuming A and B are both very optically thick. The 
model reproduces the data reasonably well at almost all 
velocities, and much more closely overall than the other 
models considered. 

The good fit of this model implies that the absorption 
in one or both of the systems is produced by many small 
subunits scattered all over the continuum source from our 
point of view. In that case, the amount of light transmit- 
ted through both systems will naturally be Ia{v) x Ib{v) 
at every velocity v (Figure [7]). Deviations will only occur 
due to statistical fiuctuations, which will be greater the 
fewer subunits there are. It is more difficult, though still 
possible, to explain the observations using two 'mono- 
lithic' systems; that is, systems in which absorption from 
the ion in question arises in a single structure along our 
line of sight spanning the range of velocities seen in the 
trough, but with physical coverage of the source which 
varies with velocity (e.g., Figure 10 of Arav et al. 1999). 
Two monolithic fiows with unblended residual intensi- 
ties Ia{v) and Ib{v) can produce any blended residual 
intensity from to min(//i(u), Ib{v)) essentially indepen- 
dently at each velocity v (Figure [7]). Thus, two mono- 
lithic fiows can explain the observations, but only if they 
just happen to overlap as a function of velocity in such 
a way as to mimic the overlap of two systems of clouds. 
Such an explanation is rather contrived, and we conclude 
instead that many small subunits exist in one or both 
absorption systems. This conclusion should of course be 
tested with observations of additional overlapping ab- 
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sorption systems in other quasars, to ensure this case is 
not a fluke. 

Note that we have not considered the effects of differ- 
ent covering factors for the continuum source and broad 
emission line region. As seen in Figure [H line emission 
is a 10% effect at best, and is not a factor at all in the 
Si IV A 1393 trough of system B. 

4.1.1. Constraints on the Outflow Subunits 

The results above suggest that the absorbers A and 
B are composed of a number of optically thick subunits. 
We now discuss what we can infer about the parameters 
of these subunits, in the limit that each subunit is so 
optically thick it can be treated as opaque. 

Assume that absorber A's residual intensity at some 
velocity, /aCv), is created by Na subunits intercepting 
our line of sight, and similarly for absorber B. When the 
two absorbers overlap along the line of sight, there will 
be TV = Na + Nb subunits along the line of sight. The 
average transmitted flux i in this case will be < i 
(1 — p)^ , where p is the average fraction of the quasar's 
emission covered by an individual subunit. 

If an average N over all velocities is well defined, the 
pixel-to-pixel variations around the average value < i > 
will be distributed with variance = + a'^ , where aj 
is the instrumental error and ct^ is given by 



^2 2 



. + (1-P) 



(4) 

For example, fixed N at all velocities would have aff = 0, 
while a Poisson distribution with an average of N would 
have (T^ = N. The intrinsic variance at fixed N and 
p, crfntrinsici causcd by the random overlap (or lack 
thereof) of N subunits of uniform projected fractional 
area a. The relation between p and a, and the form of 
crfntrinsici dcpcuds ou the shape of the subunits and of 
the quasar's emitting region. In the Appendix we give 
formulae for the cases of rectangular subunits of width 
a and unit length and of circular subunits of area a, un- 
der the approximation that the emitting region of the 
quasar is projected on the sky as a square of unit area 
and uniform surface brightness (see the discussion in the 
Appendix). In both cases, CTp oc cr^. If tJa is negligible, 
there are two unknowns (a and N) and two observables 
(< i > and a) which can be used to solve for them. 

More generally, we can constrain the subunit num- 
ber and size as follows. We have a predicted profile 
i{v) = IaIb and an observed profile /(w), both of which 
depend on velocity. In our case, the wide range of i over 
the full trough and the smooth pixcl-to-pixel distribution 
of i cannot simultaneously be reproduced at fixed N . Re- 
producing the wide range of i would require a small N , 
which would not generate as smooth a velocity profile 
as observed. Each subunit will probably have a veloc- 
ity dispersion of only ~10 kms~^ (Peterson 1997), so for 
small N strong variations in i would be seen on that ve- 
locity scale. Thus, the range in i means either N oi a 
varies with velocity, or both do. To simplify the prob- 
lem, we assume the subunits have a uniform size so that 
a is constant and ctq = 0. (This should be an adequate 
approximation if the subunits have a characteristic size 
scale.) If we then assume a value for a, we can calculate 



a predicted N for each pixel as = logi/log(l — p), 
using the expression for p{a) appropriate to the chosen 
geometry. The observed profile / differs slightly from the 
predicted profile i = IaIb, due to the intrinsic variance 
on the total covering factor of N clouds {'y'intrinsic) ^''^'^ 
to the errors on I a and Ib {(Ja and us, respectively). 
Setting CTp oc CTa = as discussed above and approximat- 
ing the variance on N as cr^ = N, we have 

2, aL„„_ + (1 -p)'^N[\nil - p)f + I^A + IWb- 

... 

The probability of observing a residual intensity I ± aj 
in a pixel, given a predicted value i and associated ai , is 



P{I±ai\i±a,) = 



1 



exp 



^2 1 



(6) 

Each pixel has a different ai which depends on the 
adopted a. To choose the best model, we find the value 
of a that maximizes the likelihood of the observations: 
L = P{Ik ± (^ik l*fc i ""ifc )• Note that a systematic er- 
ror in / (e.g., due to a continuum estimate which is too 
high or too low) will yield a systematic error in a. 

We use the velocity range -700 <v< -75 kms^i to 
calculate L, avoiding both the narrow system A' and the 
high-velocity edge of the trough from system A where 
convolution with the instrumental line spread function 
may alter the true relative absorption depths in the two 
lines of a doublet (Ganguly et al. 1999). We find a best- 
fit relative filament width w = 0.0135, with a 99.994% 
(4cr) probability range of 0.0014 <w < 0.0430. We find 
a best-fit relative cloud radius r — 0.081, with a 99.994% 
(4cr) probability range of 0.029 < r < 0.143. There is no 
statistically significant difference between the likelihood 
of the two fits. 

To convert these to physical sizes, we model the 
quasar's emission as being from a Shakura & Sunyaev 
(1973) accretion disk with viscosity parameter a = 0.1 
radiating at the Eddington limit. (We discuss the is- 
sue of coverage of the quasar's broad emission line re- 
gion at the end of the section.) For this quasar we 
estimate Mbh — 6.2 x 10^ Mq from the second mo- 
ment of its Mg II emission line and its 3000 A con- 
tinuum luminosity, using the methods of Rafiee et al. 
(2007, in preparation). For those parameters, 99% of 
the continuum emission at rest-frame 1400 A comes from 
r < IbORsch, where Rsch ^ 2GMbh/c^ = 1-8 x 10^-^ cm 
is the Schwarzschild radius of the black hole. Since the 
relative sizes derived above were referenced to a square, 
not a circle, we adopt the square that has the same 
area as a circle with radius 150Rsch, which has sides 
of length ^ = 4.8 X 10^^ cm. Thus, we find a best- 
fit filament width oi w — 6.5 x 10^"* cm, with a 4(7 
range of 6.7 x 10^^ < w < 2.1 x 10^^ cm, and a best- 
fit cloud radius r = 3.9 x 10^^ cm, with a 4(7 range of 
1.4 X 10^^ < r < 6.9 X 10^^ cm. 

These sizes, small on astronomical scales, suggest an 
origin for the subunits in the accretion disk for either 
geometry. A plausible length scale for structures origi- 
nating in an accretion disk is the scale height h of its at- 
mosphere (Equation 2.28 of Shakura & Sunyaev 1973). 

^ If the accretion disk has a strong magnetic field, the pressure 
scale height may be a less plausible characteristic length. Numer- 
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At large radii, h ~ 'dR^kTs/AGMBH'mpZo where R is 
the distance from the black hole, Tg is the disk surface 
temperature and is the disk half-thickness. (Though 
not obvious from the above, h < zq because the disk 
surface temperature is lower than its midplane temper- 
ature.) In this object, the best-fit filament width equals 
the scale height /i at r = 5500Rsch = 9.9 x 10^'' cm 
and the best-fit cloud radius equals the scale height h at 
r = 25000Rsch = 4.5 x 10^^ cm. The various parameters 
for our two geometries are summarized in Table 2. 

Strikingly, the first of those distances from the central 
source is equal to the distance the absorber must have to 
cover the emission from the quasar's broad emission line 
region (BELR). As seen in Figure 4, the line emission 
in the region of the absorption troughs reaches at most 
10% of the continuum level, and at least system A covers 
both the continuum emission region and the Si IV/0 IV] 
BELR. In other transitions, both systems at least par- 
tially cover the N V and C IV BELRs, and at least system 
A covers the Mg II BELR. Since AGN BELRs are strat- 
ified, with lower-ionization gas located farther from the 
quasar, to be conservative we assume both systems lie ex- 
terior to the Mg II BELR in SDSS J0242-h0049. We use 
a relationship between L\(3000A) and i?BELR,Mgii de- 
rived from reverberation-mapping data (Peterson et al. 
2004; Rafiee et al. 2007, in preparation) to obtain 
9.1 X IQi^ cm = 5000Rsch for SDSS 
J0242+0049. Given the ±25% la scatter in this rela- 
tionship, this distance is in excellent agreement with the 
distance required for filamentary absorber subunits to 
have widths matching the disk scale height. Of course, 
the absorber could be located at any R > -RBELR,,Mgii. 
so spherical clouds of size equal to the disk scale height 
could still match the data if the outflow arises at suffi- 
ciently large radii. 

We have outlined a consistent picture wherein systems 
A and B, whether they consist of opaque filaments or 
clouds, are launched from the accretion disk exterior to 
the Mg II BELR with a subunit size comparable to the 
scale height of the accretion disk atmosphere at that ra- 
dius. As a system accelerates, its typical density will 
decrease and its typical ionization will increase, explain- 
ing the presence of high ionization species in flows aris- 
ing from a low-ionization emission-line region. When 
the systems cross our line of sight, they have line-of- 
sight velocities of vios = —2000 km s~^ for system A and 
vios = —3600 km s"'^ for system B. For System A, \vios\ 
is comparable to the VorUtai = 2900 km s^^ expected at 
its inferred launch radius of 5500i?sc/i- For System B, 
\vios\ is larger than the Vorbitai = 1400 km s~^ expected 
at its inferred launch radius of 25000i?sc/i- The spheri- 
cal cloud dispersal time would be of order ~110 years for 
T ^ 10^ K, so the subunits will not disperse on their own 
between laimch and crossing our line of sight. However, 
partial shadowing of a subunit will produce differential 
radiative acceleration of the subunit. Substantial radia- 
tive acceleration could thus shorten the subunit lifetimes 
considerably. 

One potential complication is that the observed profile 
of the overlapping trough deviates from the multiplica- 

ical simulations of accretion disks do not yet conclusively show if 
another characteristic scale is produced by magnetohydrodynamic 
turbulence (Armitage 2004). 



tive prediction (Figure 6, bottom panel) in a manner that 
is not random on velocity scales larger than ~10 kms^^. 
However, deviations on such scales should be random if, 
as expected, the individual subunits have velocity dis- 
persions of that order. Instead, the deviations seem to 
be coherent on ~100 kms~^ scales. It may be that the 
subunits do have velocity widths of that order due to 
microturbulence (Bottorff et al. 2000). Another possi- 
ble explanation is that the outflow consists of filaments 
wherein the material is accelerated so that its line-of- 
sight velocity increases by ~100 kms~^ as it crosses the 
line of sight (e.g., Arav et al. 1999). Deviations from 
the expected profile should then persist for ^^100 kms"""^ 
instead of ~10 kms^^. As compared to a model without 
line-of-sight acceleration, there could be the same aver- 
age number of filaments, but the number would change 
more slowly with velocity (although other effects, such 
as filaments not being exactly parallel, can affect that 
as well). Observations of additional overlapping systems 
would be useful for investigating this issue. 

We note that Goodman (2003) have shown that 
thin accretion disks without winds will be unstable to 
self-gravity beyond vq^i ^ 2740{10^all/MBHf^'^Rsch 
where Ie is the Eddington ratio; using the parame- 
ters adopted herein, SDSS J0242-I-0049 has rQ=i ~ 
llOORsch- However, removal of angular momentum by 
a disk wind might help stabilize a thin disk (§4.3 of 
Goodman 2003), and there is reason to believe such a 
process operates in AGN. Reverberation mapping places 
the BELRs of many AGN at r > rQ=i, and there is ev- 
idence that BELRs are flattened (Vestergaard, Wilkes, 
& Barthel 2000; Smith et al. 2005; Aars et al. 2005) as 
expected if they are located at the bases of accretion disk 
winds (Murray et al. 1995). Furthermore, quasar spec- 
tral energy distributions are consistent with marginally 
gravitationally stable disks extending out to ~ lO^Rsch 
(Sirko & Goodman 2003). 

Lastly, we note that there is no contradiction in us- 
ing the continuum source size to derive the scale size of 
the subunits for an outflow the size of the BELR. This 
is because the continuum source has a surface bright- 
ness c:i 2100 times that of the BELR. That number is 
the ratio of the continuum flux near 1400 A in SDSS 
J0242+0049 to the Si IV/0 IV] flux, which we take to 
be ~ 9, times the ratio of the areas of the Si IV/0 IV] 
BELR and the 1400 A continuum source.^ If N subunits 
of the absorber each cover a fractional area a of the con- 
tinuum source, Nx subunits of the absorber will each 
cover a fractional area a/x of the BELR. For large N 
and small a the residual intensity of each region is equal, 
i = {1 — a)^ ~ (1 — a/x)^^ , but the variance on i from 
the BELR will be a factor ~ 0.1/.X smaller than the vari- 
ance on i from the continuum source. Thus, an absorber 
covering both the continuum source and BELR will have 
essentially the same residual intensity i and variance af 
(used to derive the absorber size constraints via Equation 
6) as an absorber covering only the continuum source. 

The size of the Si IV/O IV] BELR has been measured in only 
three AGN (Peterson et al. 2004). On average, it is comparable in 
size to the C IV BELR. We therefore use the relationship between 
Lx{l350A) and -Rbelr.civ given by Peterson et al. (2006) to 
derive iZsELR.Siiv = 4.1 x lO" cm for SDSS J0242-(-0049. 
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4.2. Possible Si IV Line-Locking 

We now return to the issue of whether systems A+A' 
and B can be line-locked. Line-locking occurs when the 
reduction in line-driving flux caused by the shadow of one 
system decelerates the other, shadowed system so that 
two systems end up with the same acceleration (which 
may be nonzero). The two systems thereafter maintain a 
constant velocity separation that keeps one system shad- 
owed (Braun & Milgrom 1989). (However, there is some 
debate in the literature as to whether line-driven winds 
are unstable to the growth of shocks (Owocki, Castor, & 
Rybicki 1988; Pereyraet al. 2004). If shocks can develop, 
they could accelerate the wind out of an otherwise stable 
line-locking configuration.) For line-locking to occur in 
an accelerating flow, there are two possibilities. System 
B could have appeared along a sightline linking the con- 
tinuum source and system A+A' at 2.0280 < z < 2.0476 
and accelerated until it reached z = 2.0280 and over- 
lapped system A-|-A' at 2; = 2.0476. Alternatively, sys- 
tem A-l-A' could have appeared at z > 2.0476 and ac- 
celerated until it reached z = 2.02476 and overlapped 
system B at z = 2.0280. 

The latter scenario can be ruled out because the great- 
est deceleration of system A+A' would have occurred be- 
fore it reached z = 2.0476, when it was shadowed by the 
deepest part of system B. Instead, the deepest part of 
system B is observed to be shadowed by the shallowest 
part of system A. If line-locking was going to occur in 
this scenario it would have had to set in when the shad- 
owing was greatest (or earlier than that, if less than full 
shadowing produced sufficient deceleration). If it did not 
happen then, it will not happen with the observed, lesser 
amount of shadowing. 

The former scenario of an accelerating system B which 
has ended up line-locked is plausible. The observed shad- 
owing as a function of velocity could in principle have 
halted system B. 

One requirement of this former scenario, however, is 
that the narrow absorption at z = 2.0476 (system A') 
should not be associated with system A, the broad ab- 
sorption immediately shortwards of it. If they were asso- 
ciated, then some of the gas in system B at —350 < z < 
—50 kms~^ should have come to a halt at kms~^, 
where the shadowing by system A' would have been 
greater than the current shadowing by system A. Sys- 
tem A' must be located farther from the quasar than 
either system A or B, in this scenario. 

The optically thickest part of system A is likely at 
—650 < V < —450 kms~^, where numerous low- 
ionization species are seen. If any gas in system B was 
observed at w < —650 kms^^, that gas would have passed 
the point of maximum shadowing without becoming line- 
locked. In fact, no gas in system B is seen a.t v < —650 
kms~^, consistent with system B being line-locked. One 
argument against this scenario is that if system B has 
been halted by the observed shadowing, gas at differ- 
ent velocities in that system has been halted by different 
amounts of shadowing. For example, gas at —200 kms~^ 
has been halted by shadowing of only ~30% of the con- 
tinuum, while gas at —450 kms~^ has been halted by 
shadowing of ~95% of the continuum. It may be more 
physically plausible to suggest that gas at —450 kms^^ 
has been halted, but that gas at —200 kms~^ has not 



yet been sufficiently shadowed to become line-locked. In 
other words, in this model system B is in the process of 
becoming line-locked. However, comparison of the SDSS 
and UVES spectra shows no evidence for variability in 
these Si IV troughs. The timescale for velocity changes 
in this scenario could be longer than 1.4 years (rest- 
frame), which would rule out line locking in a Murray 
et al. (1995) disk wind in which the entire acceleration 
phase lasts ~ 1.6 years, or the line-locking could be oc- 
curing in a helical flow, stable on timescales of years, 
in which our sightline intercepts the flow before the gas 
becomes line-locked. 

Finally, note that the Si IV profiles in SDSS 
J0242+0049 are intriguingly similar to some of the poten- 
tially line-locked N V profiles seen in RX J1230.8+0115 
(Ganguly et al. 2003). The z = 0.1058 system in 
that object has a profile similar to that of system A+A' 
(strongest absorption at both ends of the profile) , and its 
z = 0.1093 system is similar to that of system B (opti- 
cally thick, with the strongest absorption in the middle 
of the profile, at a velocity corresponding to the weak- 
est absorption in the other system). Both systems have 
only about half the velocity widths of those in SDSS 
J0242+0049, however, and the relative velocities of the 
two systems are reversed — the weaker, single-peaked 
absorption profile has the lower outflow velocity. It is 
also worth noting that the Lya absorption profile in each 
object appears to share the same covering factor as the 
species discussed above, while at least one moderately 
higher-ionization species in each object (N V here, and 
O VI in RX J1230.8+0115) has a larger covering factor 
which yields nearly black absorption troughs. Whether 
these similarities are just coincidences will require data 
on more candidate line-locking systems. (The line-locked 
systems in Q 1511+091 studied by Srianand et al. (2002) 
are much more complex, but do not seem to include any 
profiles similar to those in SDSS J0242+0049.) 

5. CONCLUSIONS 

We flnd that the C IV BAL trough at z = 1.87988 in 
the spectrum of SDSS J0242+0049 [v = -18400 kms'^ 
relative to the quasar's rest frame) has likely undergone 
an acceleration of a = 0.154 + 0.025 cm s^^ over a period 
of 1.39 rest-frame years. This is the largest acceleration 
yet reported in a BAL trough >1000 kms~^ wide. 

We also derive constraints on the outflow properties of 
two absorption systems, overlapping and possibly line- 
locked in Si IV, atz = 2.0420 and z = 2.0254 {v = -2000 
kms~^ and v = —3600 kms~^ relative to the quasar, re- 
spectively). The overlapping trough in common to both 
systems indicates that at least one of the systems must 
consist of individual subunits. This contrasts with re- 
sults strongly suggesting that the BELR itself consists of 
a smooth flow, rather than a clumped one (Laor et al. 
2006), but agrees with results for a narrow intrinsic ab- 
sorber in the gravitational lens RXS J1131— 1231 (Sluse 
et al. 2007). 

Assiiming identical, opaque subunits, our data are con- 
sistent with Hplieric;al clouds of radius r ~ 3.9 x 10"'^^ cm 
or linear filaments of width w ~ 6.5 x 10^** cm. These 
subunits must be located at or beyond the Mg II broad 
emission line region. At that distance, the above fila- 
ment width is equal to the predicted scale height of the 
outer atmosphere of a thin accretion disk. Insofar as that 
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is a natural length scale for structures originating in an 
accretion disk, these observations are evidence that the 
accretion disk is the source of the absorption systems. 
It would be useful to obtain high-resolution spectra of 
additional cases of distinct but overlapping intrinsic ab- 
sorption troughs in quasar spectra to determine if this 
case is representative. If so, it would also be worth ex- 
tending this work's analytic study of the implications of 
the residual intensity variance to numerical studies in- 
cluding a realistic quasar geometry, a range in absorber 
sizes and optical depths, etc. 
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APPENDIX 

Consider the case of an absorber consisting of opaque subunits of a uniform shape. Suppose our line of sight to a 
quasar's emitting regions is intercepted by N of these subunits, randomly distributed transverse to the line of sight. 
Then the scatter possible in the covering fraction at fixed N due to the random overlap (or lack thereof) of the subunits 
with each other will depend on the shape of the subunits. To obtain expressions for this variance, we approximate 
the quasar's emitting regions as a square of uniform surface brightness on the plane of the sky. We do this solely 
because expressions for the variance have been derived for the case of the unit square covered by two relevant subunit 
geometries: circles of area a and filaments of unit length and width a. We take the first case to represent a true cloud 
model, and the second to represent a magnetically confined 'filament' model. 

The case of the unit square randomly overlapped by filaments parallel to each other and to two sides of the square, 
and of unit length and width a, is treated by Robbins (1944). The unit square is defined as the set of points 
{0<a;<l;0<?;< 1}. The filaments that overlap the square are centered at y = 0.5 and distributed randomly in x 
over — I < X < 1 + f. Because of edge effects, the average area covered by a filament is p = and the average area 

uncovered by N filaments is i = (1 ~ p)^ . The variance in the fractional area covered is 

2a\{l - p)^+i - (1 - a)(l - 2p)^+i] 



Ornament. = (1 - - 2p)^ - (1 - pY" + 



\2N 



{N+l)p 

2a^[il - p)^+^ - jl - 2p)^^ 
(7V + l)(iV + 2)p2 



(1) 



for a < 0.5. 

In the case of the unit square randomly overlapped by circles of area a, circles that overlap the square are distributed 
such that their centers are within a distance r = y^a/n of the unit square. Again the average area uncovered by N 
circles is given by i = (1 — p)^ , but in this case p = 7rr^/ (1 + 4r -I- 7rr^). The variance in the fractional area covered 
can be derived from expressions given by Kendall & Moran (1963), yielding 



^circles 



1 -I- 4r 



1 -I- 4r 



, 64 
1 - Anr"' H r 



8r^ I - 



1 -t- 4r 



1 + 4r 



2N 



2r 



1 - 



2r^(7r — cos 



2r 



^ sin (cos 1 ^) 



2r 



1 + 4r 



N 



{-nq - Aq^ + q^) dq 



(2) 



for a < 0.5. The integral must be evaluated numerically for most N. 

For the same a and N, cr'^^^^^^s > '''filaments- This can be understood by placing a subunit of either type in the 
center of the square and considering the probability that a second subunit of the same type will overlap the first. 
There is an area 2a in which a second filament can be placed to have some overlap with the first (filament centers at 
0.5 — a < X < 0.5 -|- a). There is an area 4a in which a second circle can be placed to have some overlap with the first 
(circles centered within of {0.5,0.5}, for an area of 7r(2-y/^)^ — 4a). If a is small, the most likely value of i is 

i = \ — 2a for both geometries, but with circles there is a higher probability of i > 1 — 2a and thus a larger variance. 
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TABLE 1 

SDSS J0242+0049 Spectroscopic Observations and Inferences 





SDSS 


SDSS 


Epoch 


^trest 


Si IV, C IV Shift 


Si IV, C IV Shift 


Source 


Plate 


Fiber 


in MJD 


(days) 


vs. MJD 52188 


vs. MJD 53619 


SDSS (1) 


408 


576 


51821 


-80 


0, 


1,4 


SDSS (2) 


707 


332 


52177 


36 






SDSS (3) 


706 


617 


52199 


43 






SDSS Avg. (2+3) 






(52188) 


40 




1, 3 


SDSS Avg. (1+2+3) 






(52066) 







1, 3 


UVES 






53619 


507 


1, 3 





Note. — Epochs are given on the Modified Julian Day (MJD) system. The rest-frame time interval Atrest is given relative to 
MJD 52066. Velocity shifts of absorption lines are given in SDSS pixels (69 kms~^); the C IV shift is the first number and the 
Si IV shift is the second number. 



TABLE 2 

SDSS J0242+0049 Sdbunit Parameters 



Subunit 
Geometry 


Avg. Number 
of Subunits JV 


Best-fit Relative 
Width or Radius 


Relative 99.994% 
Confidence range 


Best-fit Physical 
Width or Radius 


Physical 99.994% 
Confidence range (cm) 


Atmospheric 
Scale Height Distance 


Filaments 
Spheres 


203 ± 81 
177 ± 71 


0.0135 
0.081 


0.0014 - 0.0430 
0.029 - 0.143 


6.5 X 10" cm 
3.9 X 10^5 cm 


6.7 X 10^3 - 2.1 X 10'^'' 
1.4 X 10^5 - 6.9 X lO^"' 


9.9 X 10^'^ cm = 5500 Rs^h 
4.5 X 10^** cm = 25000 Rsch 



Note. — The average number of subunits N is the number of subunits responsible for absorption at each pixel, averaged over all pixels. The total 
number of subunits present depends on the unknown velocity width of each subunit. The atmospheric scale height distance is the distance from the 
black hole at which the accretion disk atmospheric scale height equals the best-fit width or radius of the subunit in question; see §4.1. Rsch refers to 
the Schwarzschild radius of a black hole with mass 6.2 x 10® M©. 



4000 6000 8000 
— \ ^ \ ^ \ ^ \ ^ \ — ^ — ' 




1500 2000 2500 3000 

Wavelength in A (log scale) - Top: Observed - Bottom: Rest at z = 2.062 

Fig. 1.— VLT UT2 + UVES spectrum of SDSS J0242+0049, smoothed by a 1 A boxcar filter. 
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2000 
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^ -4000 -2000 2000 4000 

Broad Absorption from Lya and NV, SilV and CIV Doublets vs. velocity at z= 1.87988 

Fig. 2.— UVES spectra of BAL troughs in SDSS J0242+0049 vs. velocity (in kms^^) in the z = 1.87988 frame. Negative velocities 
indicate blueshifts and positive velocities indicate redshifts relative to that frame. Zero velocity corresponds to the long-wavelength members 
of doublets, and dashed vertical lines indicate all components of each transition. Contaminating narrow absorption lines are present near 
all troughs, but especially in those found shortward of the Lya forest. 
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CIV shifted 



0„ 



-6000 -4000 -2000 2000 4000 6000 

Average SDSS and UVES spectra for z= 1.87988 CIV and SilV BALs vs. velocity (km/s) 



Fig. 3. Comparison of the z = 1.87988 C IV BAL in SDSS J0242+0049 at the average SDSS epoch and the UVES epoch. Negative 
velocities indicate blueshifts and positive velocities redshifts, relative to z = 1.87988. The solid line is a weighted average of all three SDSS 
spectra. The dashed line is the UVES spectrum binned into the same pixels as the SDSS spectra. Dotted vertical lines indicate the fitting 
regions used when conducting the test. The top panel compares the unshifted spectra for the Si IV trough, and the middle panel the 
unshifted spectra for the C IV trough. The bottom panel compares the C IV troughs after shifting the average SDSS spectrum toward 
shorter wavelengths by 3 pixels. 
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1380 1385 
H — I — I — I — I — I — I — ^ 



1390 1395 
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1400 
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1410 




Observed 
<- Absorption 
Redshifts 



Si IV Emission at z=2.062 
<- Expected Line-Locking Redshifts 



4200 4210 4220 4230 4240 4250 4260 4270 4280 4290 4300 
SilV Doublets in SDSS J0242+0049 - Bottom: Observed X (A) - Top: Rest X (A) at z=2.0476 



Fig. 4. — Two broad, overlapping Si IV doublets in the unnormalized spectrum of SDSS J0242+0049. Line identifications and redshifts 
for the different troughs are given on the figure. There is also narrow Si IV absorption at z=2.0314 which is not marked. 
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SilV in SDSS J0242+0049 - Normalized Flux vs. Av (km/s) from z=2.0476 or z=2.0280 

Fig. 5. — Velocity plot of Si IV absorption after normalization by a fit to the total spectrum (continuum and weaJf emission lines). 



Quasar Outflow Constraints 13 




Velocity (km s"') 



Fig. 6. — Fits to the blended Si IV trough. The trough containing blended absorption from both redshift systems is shown as the solid 
line in all panels. The fits are shown as lighter lines with total error bars that include the observed errors on the flux in the blended trough, 
so that at each pixel the deviation between the actual trough and the fit can be directly compared to the total accompanying uncertainty. 
Top panel: all three observed Si IV troughs are ovcrplotted. The dashed line shows the unblended trough A1393 trough, plotted in the 
z = 2.0280 frame. The dot-dashed line shows the unblended trough A1402 trough, plotted in the z = 2.0476 frame. Second panel: the 
fit and errors shown arc for an optically thin lowcr-rcdshift system. Third panel: the fit and errors shown arc for an optically thick 
lowcr-rcdshift system with maximum overlap in covering factor with the optically thick highcr-rcdshift system. Bottom panel: the fit 
and errors shown arc for the case where each system's covering fraction describes its fractional coverage of the other absorption system, so 
that the residual flux from both optically thick troughs can be multiplied together to give the residual flux in this blended trough. 
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I,(X) = 0.5 Ib(X) = 0.5 



g 1,JX) ^ 0.5 





<Ia.b(^)> 

= <IaXIb> 
= 0.25 



Fig. 7. — An illustration of how observations of optically thick absorption systems overlapping in velocity space can constrain absorber 
substructure. Each square is a schematic depiction of a quasar's emission region at the same wavelength A. The hatched zones are the areas 
of the emission region covered by an absorber. The leftmost colunm depicts absorption by system A only, which is assumed to produce a 
normalized residual intensity of Ia(-^) = 0.5 at the wavelength shown. The value of Ia(-^) is the same regardless of whether the absorption 
is due to a monolithic flow (top row) or due to randomly placed subunits (bottom row), here shown as spheres in projection. Similarly, 
the second column from the left depicts absorption by system B only, which is also assumed to produce a normalized residual intensity of 
Ib{^) = 0-5 at the wavelength shown. When systems A and B overlap, monolithic flows can produce any normalized residual intensity in 
the range < Ia+b{^) < 0.5 (top row, right) depending on the specific areas covered by each flow at that wavelength A. On the other 
hand, flows consisting of many randomly placed subunits will naturally yield an average value of </a+s(-^)>=<^a(-^) x Ib{-^)>= 0.25 at 
every wavelength A where the systems overlap (bottom row, right). Deviations from this average will occur due to statistical fluctuations, 
which will be smaller the more subunits there are. Measuring the deviations thus enables us to constrain the size and number of the 
subunits (see §4.1.1). 



